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Towards a strong-coupling theory of QCD at finite density* 

B. Bringoltz and B. Svetitsky^ 

'^School of Physics and Astronomy, Raymond and Beverly Sackler Faculty of Exact Sciences, Tel Aviv 
University, 69978 Tel Aviv, Israel 

We apply strong-coupling perturbation theory to the QCD lattice Hamiltonian. We begin with naive, nearest- 
neighbor fermions and subsequently break the doubling symmetry with next-nearest-neighbor terms. The effective 
Hamiltonian is that of an antiferromagnet with an added kinetic term for baryonic "impurities," reminiscent of 
the t-J model of high-Tc superconductivity. As a first step, we fix the locations of the baryons and make them 
static. Following analyses of the t-J model, we apply large- methods to obtain a phase diagram in the Nc,Nf 
plane at zero temperature and baryon density. Next we study a simplified 17(3) toy model, in which we add 
baryons to the vacuum. We use a coherent state formalism to write a path integral which we analyze with mean 
field theory, obtaining a phase diagram in the (ns,T) plane. 
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Color superconductivity |^ at high density is 
so far a prediction only of weak-coupling analy- 
sis, valid (if at all) only at very high densities. We 
seek confirmation from methods that do not de- 
pend on weak coupling, as well as an extension to 
the regime of moderate densities. Since Euclidean 
Monte Carlo methods are unavailable when the 
chemical potential is non-zero, we turn to the 
strong-coupling limit of QCD; harking back to the 
early days of lattice gauge theory, we work in the 
Hamiltonian formulation. We derive an effective 
Hamiltonian for color-singlet states that takes the 
form of an antiferromagnet with a kinetic term for 
baryons. This effective Hamiltonian is very diffi- 
cult to study. As a first step, we fix the position 
of the baryons and study mesonic excitations in 
the baryonic background. Coherent-state meth- 
ods then enable us to derive equivalent models 
that are tractable in various limits of large Nc 
and/or Nf. We benefit from the considerable 
work done on antiferromagnets in the context of 
large-Tc superconductors as well as from 

the early work of the SLAC group |] and of Smit 
|6j on strong-couphng QCD. 
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1. The effective Hamiltonian 

The lattice gauge Hamiltonian is composed of 
electric, magnetic, and fermion terms, 



H — He + Hjj + Hp , 



(1) 



where the first term is the unperturbed Hamilto- 
nian in the strong coupling limit. 



(2) 



The ground state sector of He is highly degen- 
erate, consisting of of all states with zero electric 
flux, whatever their fermion content, 
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(3) 



Neglecting the magnetic term, which only con- 
tributes in high order, we perturb with the 
fermion kinetic term. 



H. 



j>0 



(4) 



We use four-component fermions with a general 
(diagonal) kernel D{j). This is the best we can 
do, since domain- wall fermions are unavailable to 



2 



us at strong coupling [Q and there is no Hamil- 
tonian overlap formalism |^ . We will discuss the 
properties of these fermions in a moment. 

The degeneracy is lifted in second order via di- 
agonalization of an effective Hamiltonian, 

x(^'^^'V)„+,.,, (5) 

where sj^ is a sign factor and K{j) is a new ker- 
nel. The effective Hamiltonian contains fermion 
bilinears at each site, 

Q;1= (V'^MV)„, (6) 

where the Af matrices act on the Dirac and the 
flavor indices, = (g) A". The operators Ql 
generate a U {4:Nf) algebra. If we fix the baryon 
number i? on a site, the color-singlet states on 
that site make up an irreducible representation 
of this algebra, whose Young tableau has m = 
B + 2Nf rows and Nc columns (see Fig. 1). 
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Figure 1. Young tableau for single-site states with 
B = m-2Nf. 



Expressed in terms of Q^, the effective Hamil- 
tonian is an antiferromagnet. If the kernel D(j) 
is chosen to give nearest-neighbor couplings only, 

(2) 

-ff^ff will take the form 

^^iff^=^(l)EQM+A- (7) 

This antiferromagnet has the accidental U{ANf) 
symmetry of naive fermions, which is responsi- 
ble for part of the doubling problem. The even-j 



terms in (^) break this symmetry to SU{Nf) x 
SU{Nf) X U{l)v X U{1)a, which [apart from the 
unbreakable axial J7(l)] is the desired continuum 
symmetry. Since we are interested in strong cou- 
pling where the free fermion dispersion relation is 
of no interest, this might be a good-enough par- 
tial solution of the doubling problem. 

(2) 

The theory of H^^f' contains only static 
baryons, which make their presence felt through 
fixing the rep of U{ANf) on each site. These 
baryons can move in the next order in perturba- 
tion theory (only if Nc = 3 — a fortunate special 
case) . The effective Hamiltonian in third order is 

+/I.C., (8) 

where the baryon operators are the color singlets 
— ea/37V'S^^V'7 that belong to the II rep 
of U{ANf). The simple form of iJofr is deceptive, 
for the 6^'s are composite and hence do not obey 
canonical anticommutation relations. 

For nearest-neighbor fermions, the usual spin 
diagonalization |^ gives a simplified Hamiltonian, 

The complete Hamiltonian H]^g + H^g resem- 
bles that of the t-J model, which represents the 
strong-binding limit of the Hubbard model and is 
much studied in connection with high- Tc super- 
conductivity But our Hamiltonian is much 
more complex. 

2. Static baryons 

Let us beat a strategic retreat to the second- 
order theory, where baryons constitute a static 
background. If we begin with the nearest- 
neighbor model, in a state with no baryons, then 
the effective Hamiltonian (^ is that of a U{4Nf) 
antiferromagnet with spins in a rep specified by 
Nc and by m = 2Nf (which can vary from site to 
site). This can be studied in the limits of large 
Nf , Nc by various transformations , and 

the result — still for B — — is the phase diagram 
in Fig. 2. (Note that the coupling constant K{1) 
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Figure 2. Phase diagram for the nn antiferromag 
net at T = 



is just a scale, which has no effect at T = 0.) 
The large-iV/ phase is disordered, and wiU not 
concern us. The location of the phase boundary 
can be established by studying a Schwinger boson 
representation of the U{4Nf) spins, and its slope 
turns out to be N^/Nf = 0.31; this means that 
the QCD vacuum is safely in the ordered phase 
for any reasonable number of flavors. 

The ordered phase is conveniently studied in a 
cr model representation, which comes from rewrit- 
ing (|^) in a basis of spin coherent states |^ . (This 
is valid for any iVc, Nf but proves soluble in the 
Nc oo limit.) The degrees of freedom of the a 
model are the N x N matrices (with N = 47V/) 



where 



A: 



l-Ar/2 I 

-1 



N/2 



(10) 



(11) 



The field Un runs over the group U{N), and 
the manifold covered by Qn is the coset space 
UiN)/[U{N/2) X U{N/2)]. The action of the a 
model (in continuous time) is 



dT 



where, in terms of the matrices Qn, the nearest- 
neighbor Hamiltonian takes the form 

H=['-^\ A'(l)y"TrQnQ„+^. (13) 



Clearly as iVc ^ oo the ground state is the 
classical minimum of iJ, in which on all sites 
Qa = Qo, which can be rotated to Qn = A. 
Thus the symmetry is spontaneously broken as 
U{N) U(NJ2) X U{N/2), with N^/2 Gold- 
stone bosons g. 

We now restore next-nearest-neighbor cou- 
plings, viz. 



H' 



K{2) 



(14) 



where Ql = 2Tr Af''Q„. The symmetry of the 
theory, as discussed above, is U {Nf) x U (Nf); the 
classical minimum is at Qn = 70 , which breaks all 
the axial generators and leaves the vector U{Nf) 
unbroken. This is what we would expect for the 
ground state in the vacuum sector. 

3. Adding baryons 

The B = states considered above were spec- 
ified by choosing the m — 2Nf rep on each 
site. Choosing a different rep adds (or subtracts) 
baryons on a site-by-site basis. For instance, we 
can add a single baryon by adding a row to the 
Young tableau (Fig. 3). One can similarly add 



(12) 



B=0 B=0 B=l 8=0 B=0 



Figure 3. Adding one baryon to a site in the 
Nf = 1 theory [a U{4) antiferromagnet] 



baryons on an entire sublattice. The TVc — > 00 
limit directs us to find the classical ground state 
of the Hamiltonian, which always breaks the sym- 
metry spontaneously along the lines shown in the 
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section above. We can study the effects of the 
fluctuations by doing mean field theory for finite 
Nc- To do this, we drop the kinetic term in ([l^ ) 
and go to T ^ by calculating the resulting clas- 
sical partition function. 

4. Mean field theory 

In MF theory, we write down a trial Hamilto- 
nian and calculate a variational free energy $. 
The simplest trial Hamiltonian, containing no 
site-site correlations, is 



Hq — Qn ■ hn , 



(15) 



where the magnetic fields are variational pa- 
rameters. (We write Qn for the vector whose 
components are Q^.) The free energy obeys 
F < $ with 



$ = - XI 



(mn) 



(16) 



Here K = P^^-^ i'^Y ^^'^ magnetiza- 



tion. 



^0„ 



(17) 

Note that the integration measure dQn depends 
on the U{N) representation chosen for site n. 

We minimize $ with respect to {/in} and get 
the mean field equations, 

. JdQne-"^-^--^ 



M: 



'Qn 



(18) 



/dQne-^'5"-S"^(")'''" 

If all sites are in the same rep of U{N), then the 
bipartite nature of the antiferromagnetic system 
gives two coupled sets of MF equations. Other- 
wise one gets as many coupled MF equations as 
there are inequivalent sites. 

After solving for fl^ we evaluate $ via 



Elog i^j dQe 

-K E Mn • Mn 

(mn) 



(19) 



We seek the global minimum of $. Once we find 
it, we can examine its symmetry properties and 
identify the phase favored at temperature T . 

5. ?7(3), a toy model 

The Nf = 1 theory, the U{A) antiferromagnet 
(with impurities), contains many degrees of free- 
dom in which to do mean field theory. A simpler 
non-trivial model is a toy model with U{3) sym- 
metry, which does not correspond to an actual 
value of Nf. The most symmetric B — state 
in this model cannot have i? = on every site, 
but must alternate between the reps correspond- 
ing to B = ±1/2 (see Fig. 4). A B ^ state 



B=-l B=l B=-l 8=1 B=-l 



Figure 4. The B = state in the f/(3) toy model 



is specified by breaking the alternating pattern, 
as in Fig. 5. We in fact create a non-zero den- 



3=-l B=l B=+l B=l B=-l 



Figure 5. Adding a baryon in the /7(3) toy model 



sity of baryons by adding a row to one or more 
sublattices, forming a lattice with a unit cell of 
2x2x2 sites. The unit cell may contain sites 
with TO = 1, 2, or even 3. For to = 1, 2 the mani- 
fold of Qn is C/(3)/[C/(2) X [/(I)], while for to = 3, 
the singlet state, we have Qn = 0. 
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The form of [/„ for a site with m 2 is [/„ 



e with 









0-1 








02 


-al 


-0-2 






(20) 



We write 

ai = pcos(6l/2)e*"i 
02 = psm{e/2)e'"\ 



(21) 



with < p, 61 < TT and < ai,2 < 27r. With these 
definitions, Q is given by 



= -2 



1 

'76 



sm p sm y 
— sin^ /9 sin 
sin^ p cos 



COS(p 

9sin( 



V 



sin2/cicos0/2 cos + </*)) 
- sin 2/9 cos 6*72 sin ( i (V' + (^)) 
sin2psin0/2 cos (^("0 — 0)) 
■ sin 2p sme/2 sin (i (ip - (j))) 



cos^ p - 



(22) 



Here </> = ai — ck2 and ip — ai + a2, with < '0 < 
47r and < < 27r. The induced measure on the 
four-dimensional manifold turns out to be 



dQ = (1 — cos^ p) cos p d{cos p)d{cos 9) 



dip dip 



(23) 



27r 4tt 

An example: For the B = case shown in 
Fig. 4, the inequivalent sites are just the even and 
odd sublattices. The corresponding MF equa- 
tions are 



Ml 



M2 



/ dQ e+^^'^ P^Q 
JdQe+<^^QP^ 



(24) 



Cases with B ^ will give more coupled sets 
of MF equations, according to the number of in- 
equivalent sites in the unit cell. 

For each baryonic configuration we obtain a 
phase transition as a function of temperature. In 
some cases it is first order and in others, second 
order. In all cases the symmetry breakdown at 
low temperature is U{3) ~> U{3)/[U{2) x C/(l)]. 
As we increase the baryon density the transition 



temperature decreases, but it never vanishes; MF 
theory always breaks symmetry at T = 0, even in 
one dimension. We summarize our findings in the 
phase diagram in the temperature-density plane 
shown in Fig. 6. 

The future holds, we hope, the removal of the 
various approximations that led from QCD to the 
U (3) toy model. To begin with, we must do better 
than classical mean field theory and include the 
quantum kinetic term in the a model (p^. The 
U{3) model must be generalized to the J7(4iV/) 
symmetry group of naive fermions, which should 
be broken to U{Nf) x U{Nf) by the nnn cou- 
pling. Going beyond the static baryon picture, we 
can disorder the baryon background by a replica 
method; eventually dynamical baryons should be 
included with the third-order kinetic term. Once 
the theory becomes realistic enough, we can com- 
pare the results at each step to the weak-coupling 
predictions of color superconductivity with vari- 
ous values of Nf [|lo| . 

There is, however, a limitation inherent in the 
strong-coupling theory. As we saw in the mean- 
field analysis above, one is easily misled into fil- 
ing the lattice with baryons to saturation. The 
baryon density is limited according to 



(25) 



B 2Nf 
— < 

and strong coupling means large lattice spacing 
a. It is possible that color superconductivity will 
not show up at any density short of saturation. In 
that case we will have to content ourselves with 
a new effective theory for baryonic matter, short 
of the transition to quark matter. 
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Figure 6. Classical phase transition of the U{3) antiferromagnet with a density of baryon impurities. The 
ordinate is proportional to temperature, while the abscissa is baryon density in units of the maximum 
saturation density Pmax- 
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